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Abstract 

The SLR technology experienced a large number of technical improvements 
since the early 1970 les , leading now to a millimetric instrumental accuracy. 
Presently it appears as useless to increase these instrumental performances as long 
as the atmospheric propagation delay suffers its actual imprecision. It has been 
proposed since many years to work in multiwavelength mode, but up to now the 
considerable technological difficulties of subpicosecond timing have seriously 
delayed such an approach. 

Then a new possibility is proposed, using a device which is not optimized 
now for SLR but has already given good results in the lower troposphere for wind 
measurement: the association of a radar and a sodar. While waiting for the 2-X 
methodology, this one could provide an atmospheric propagation delay at the 
millimetre level during a few years with only little technological investment. 


1/ INTRODUCTION 

It has been pointed out since a long time that all space geodesy techniques 
have to deal with the same general problem, i. e. the crossing of the atmosphere to 
reach either an artificial satellite, or the moon, or a star, or a quasar, etc... 

Nevertheless, if it is the same atmosphere for every techniques, the effect is 
known to be quite different for radiowave and for optical methodologies, the first 
ones suffering more than the other ones from the crossing of ionosphere (but this 
is corrected classically by two-frequencies methods) and from the troposphere 
transit (because of the atmospheric water vapor content which is quite 
unpredictable and has a strong impact on propagation of radio waves and not on I 

optical ones). Considering these aspects, the general advantages of radio techniques I 

compared with optical ones in geodesy are their all-weather capabilities, and their J 

drawbacks are linked with their poor tropospheric correction quality. 

As long as the instrumental accuracies were at a few centimeters level, this 
meteorological aspects were of secondary concern. This is no longer the case since 
many years with VLBI, who had to use costly water vapor radiometers to upgrade 
the tropospheric correction to an excellent level (and since this period, VLBI has got _ 

the best positioning precision among all space geodesy techniques). Since a few 
years too it is no longer the case with satellite laser ranging, which is able of an 
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internal instrumental precision of a few millimeters, to be compared with the final 
positioning accuracy of 2-3 cm. Of course the tropospheric corrections are not the 
only problems biasing the results (consider for example the poor world coverage of 
SLR stations), but significant improvements of this parameter must be researched. 
Many teams have pointed out that multiwavelength methodology could pnivide this 
required amelioration. Anyway it is quite clear that this will ask for a tremendous 
technological effort and will probably not be operational since some years. For that 
reason we have looked for a new solution, able to give better tropospheric correction 
with up-to-date techniques, even if probably not as efficient as 2-colour method, 
but immediately available. 


n/ THE ATMOSPHERIC CORRECTION FOR SLR 

If n is the index of refraction, n = c 0 / c (c 0 is the light celerity in vacuum 
and c in the atmosphere). 

Classically we define the co-index of refraction as 

N = ( n - 1 ).10 6 

and for optical wavelengths, using for example Essen's formula: 


N = A, ( P. / T ).{ 1 + ( a - b.t ).P. } - B x ( P v / T ) 


where P, is the atmospheric pressure, P v is the water vapor pressure. T is the 
absolute temperature and t is the centigrade temperature. A x and B x are \- 
dependant parameters, and a and b are constants. With a 10' 7 precision it is 
acceptable to use the simplified formulation: 

N = A\ ( P, / T ) - B x ( P v / T ) 


We define also the geometric length of the optical ray L g , the optical path L„ , 
and the geometric distance L between two points s 0 and s, being one at the ground 
level and the other beyond the atmosphere, at the satellite level: 


L g = fds L. = f n(s).ds 

■*6 *0 

If we call AL = L 0 - L , it is a function of the angle a between the vertical and 
the direction of the satellite. We notice that: 


AL = ( L 0 - L g ) + ( L g - L ) 


As the curvature of the ray path is low. and as this curvature is experienced 
only on a small range (a few tens of kilometers), the difference L ( - L is generally 
considered as negligible. On another hand. AL can be expressed as: 

AL = f(a) . AL vertlcal 
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And the function f(a) may be found in (Berrada-Baby et al., 1987 or 
Akhundov et Stoskli, 1992). We have now to deal with AL for 4 = 0 . ho being the 
elevation of the SLR station above "sea level", 


A L = [ N(h).dh 
*0 

Considering the very low dependance of the result regarding the water vapor 
for these optical wavelengths, we will focus on the "dry" part of this expression AL d : 


M d 



.dh 


And, p(h) being the air density, T(h) the absolute temperature. P(h) the 
atmospheric pressure and g(h) the acceleration of the pesanteur at the elevation 
h, with R = 287 J.kg/°K 

* P. = p(h) . R . T(h) 

* d P(h) = p(h) . g(h) . dh 

so that. 


A L d 



dPjh) 

8 (h) 


.A . R 


At this level we may have different assumptions: 

1 st hyp. g(h) = go = constant. Then AL d = A.R.P 0 / & 

2 nd hyp. g(h) = go / ( 1 + 2.h/r 0 ) closer to the reality. We find then: 


M d = 


A.R.Pn 


So 


1 + 


2/?0 


'o-'ir ( 80 + ^*P) - 2 ,g 0 


( p is the gradient of T ). For P 0 = 1000 mbar, we find AL d * 2,3 meters 


And the difference between these two models is below 1 cm : various studies 
show (Berrada-Baby. 1987) that this discrepancy is around 5 mm. and Is quite 
stable and thus easy to model. On another hand, the differences between the 
second model and real data deduced from (Cira 1965) campaigns are quite small 
(1 mm typically). 
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So we conclude that it is extremely important to measure P 0 with an 
excellent precision (1 mbar of error induces 2.3 mm on the atmospheric correction). 
And one must take into account the fact that the function ffa) (that is the air-mass 
number) will multiply these values by numbers up to 2.5 in operationnal SLR 
measurements . 

If the temperature is measured in the first 5 kilometers, looking at fig. 1 and 
fig. 2, we observe that half the total correction Is already acquired (fig. 1), and that 
the upper layers are quite well defined by the profile In the troposphere. If we 
compare these facts with the abovementionned value of 5 mm, whose noise is 
around 1-2 mm, multiplied by the air-mass value (1 to 2.5), i. e. less than 5 mm, 
it is clear that any temperature profile of good precision acquired in the troposphere 
will leave a residual error on the atmospheric correction at the millimeter level. 



Fig. 1 : Fraction (Y-axis) of the total atmospheric correction for a vertical transit 
of a laser pulse, from the ground level up to a given altitude (X-axis, in 
meters). Computed from an observed radio-profile in Greece, 1990. 
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in/ PROPOSITIONS 


The first one Is obviously to measure P 0 quite carefully. It is necessary to 
calibrate often the barometers employed, to measure preferably the pressure close 
to the level of the axes intersection of the telescope (it is generally the case). It must 
be possible to measure the pressure with an absolute precision < 0.2 mbar. 


Considering the efficiency of barometric levelling, whose precision may reach 
one meter , we observe that it means that constant pressure surfaces are quite 
horizontal, so that it is useless to measure the horizontal gradient of P 0 in order to 
correct its effect in the direction of the sight, provided the weather is reasonably 
quiet. 

Anyway, it seems mostly advisable to improve the correction by a good 
measurement of T(h) in the direction of the satellite, at the 0.1°C level, up to an 
elevation' of 3 to 10 km. 

In these conditions, one may be sure to get an atmospheric correction better 
than 1 mm. It is easy to notice that such a precision with 2-wavelength SLR will 
require a 20 times better precision on the differential time-of-ffight measurement 
between the two colors (i. e. 0.5 picosecond), which is quite uneasy to reach. 

The solution we propose for an easy measurement of the temperature profile 
along the line of sight is a SODAR (Acoustic LIDAR) associated with a RADAR, the 
acoustic and the radio wave having the same wavelength. The radar is used to track 
the acoustic wave, so it allows to measure the speed of the sound in the 
atmosphere, which in turn provides an excellent temperature profile (at the 0.1 °C 
level). Such instruments (e. g. Remtech) are now used close to airports, in order to 
measure winds and wind shears, and this technology is already available. The 
common ranges are 3 to 5 km. but large instruments may reach 10 km. 
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TV/ CONCLUSIONS 

We propose that an experiment be carried with a radar/ sodar equipment 
close to an^ up-to-date SLR station in order to check the possible 1 mp r <^emen 
that such data could provide to the laser data. This solution is probably 
^cellent alternative allowing to wait for the 2-wavelength generation of SLR 

stations. , 

On another hand, we insist on the necessity to check quite c^ef^tythe 
barometric equipments used in SLR stations, as the P. parameter is by far the most 
important to perform a good atmospheric correction. 
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